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Abstract 

The dimuon invariant mass spectrum is searched in the range between 5.5 and 14 GeV 
for a light pseudoscalar Higgs boson a, predicted in a number of new physics mod- 
els, including the next-to-minimal supersymmetric standard model. The data sample 
used in the search corresponds to an integrated luminosity of 1.3 fb -1 collected in pp 
collisions at a/ s = 7 TeV with the CMS detector at the LHC No excess is observed 
above the background predictions and upper limits are set on the cross section times 
branching fraction a ■ B(pp — > a — > u + }i~) in the range of 1.5-7.5 pb. These results 
improve on existing bounds on the abb coupling for m a < /Wy(iS) an d are the first sig- 
nificant limits for m a > m Y (3s)- Constraints on the supersymmetric parameter space 
are presented in the context of the next-to-minimal model. 
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Low energy supersymmetry (SUSY) is an elegant solution to the hierarchy problem that arises 
in the Standard Model (SM), provides a candidate for dark matter, and allows for the unifica- 
tion of gauge couplings at the grand-unified-theory (GUT) scale [1-5]. However, the minimal 
supersymmetric model (MSSM) has an ad hoc Higgs superfield mixing parameter (ji) and re- 
quires very large masses for the supersymmetric partner of the top quark (stop) in order for 
the lightest CP-even Higgs boson to be heavier than 122 GeV without large stop mixing (6|. 
Both problems are solved in the next-to MSSM (NMSSM) (a review can be found in Ref. |7|), 
which extends the MSSM by introducing a complex singlet superfield which necessarily con- 
tains a scalar field component. Associated super- and scalar-potential terms generate an effec- 
tive }i parameter and easily raise the mass of the light Higgs boson without requiring a heavy 
stop 18. 9J. The added scalar field expands the Higgs sector to three CP-even scalars (hi, h2, 113), 
two CP-odd scalars (a\, aj) and two charged scalars (H + , H ). The a\ is a superposition of the 
MSSM doublet pseudoscalar (amssm) and the additional singlet pseudoscalar of the NMSSM 
(as): «i = cos 6*A fl MSSM + sin 6a<^s, where 9 a is the mixing angle. The NMSSM has two symme- 
tries that, if imposed (e.g. at the GUT scale), imply that small m ai , even m ai < (where is 
the B meson mass), and | cos 6>a| <C 1 are very natural possibilities IITOl . However, the reduced 
couplings C Hib j- = C ai}l + }1 - = C aiT + T - = tan f> cos 6a can be sizeable for large values of tan /3, the 
ratio of neutral Higgs field vacuum expectation values, even if cos 9a is small. More generally, 
superstring modeling suggests the possibility of many light a particles, at least some of which 
couple to , t + t~ and bb (ill. h" 1 the following, a (a\) denotes a general (NMSSM) light 
pseudoscalar Higgs boson. 

Searches for a light a are mainly sensitive to C flbb ffl2l [13 | . For m a < m Y (is)/ the strongest 
constraints on C flbb are those from BaBar (141H51 . For m a > mypg), only the Tevatron and Large 
Hadron Collider (LHC) have sensitivity (1611 , using production via gg — > a, where the coupling 
C fl gg derives from quark (especially bottom and top) triangle loops. This process, plus higher- 
order corrections, leads to a large cross section due to the large gg parton luminosity at small 
gluon momentum fractions, provided the C flq q (q = t,b in particular) couplings are not too 
suppressed. This large cross section will typically lead to a significant number of gg — > a — > 
events even though B(a — > }i + f.i~) is small. 

In the NMSSM context, where C flibb = tan /3 cos 9a, the existing limits HZHIll translate to rather 
modest limits on | cos 9a | • Such bounds do not strongly constrain NMSSM models of interest 
for possibly hiding a light Higgs boson because of h — > aa decays (with a — > It, 2g, 2c, 2s 
decays being dominant [19]) that are not excluded by large electron-positron (LEP) collider 
experiments |[20ll2T| . 

At tree level, the branching fraction for a — > fi + }i~ depends on m a and on tan/5, but not on 
cos $a[16]. It is nearly constant for m a > 5 GeV and ranges from 10~ 3 to 4 x 10~ 3 for tan f> = 1 
to tan j6 = 50, changing very little once tan j5 > 2. In contrast, c(gg — y a) increases rapidly with 
tan /3 due to the fact that C flbb <x tan /3. However, top-quark loop contributions and higher-order 
corrections imply a slower c(gg — > a) increase than tan 2 /3. In the context of the NMSSM, all qq 
couplings of the a\ are proportional to cos 9a, implying that u(gg — > a.\) oc cos 2 0a- 

This Letter presents the results of a search in pp collisions for a light a with a mass near the 
Y resonances decaying into two oppositely charged muons. Data used for this analysis were 
recorded by the Compact Muon Solenoid (CMS) detector in pp collisions at a center-of-mass 
energy of 7TeV, between August and November 2011. The sample corresponds to a total inte- 
grated luminosity of 1.3fb~ , collected with a dedicated trigger. As estimated in Ref. (16| and 
explicitly demonstrated here, CMS has sensitivity beyond the BaBar and CDF limits, for the 
latter due to the higher production yield [c^hcIpP ~~ ^ a ) ~ 4.5t7 Tevatron (pp — y a)] and the higher 
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acceptance and efficiency of the muon detector. Furthermore, the CMS analysis can extend the 
limits into the m a > wz Y (3s) mass range. 

The central feature of the CMS detector is a superconducting solenoid, of 6 m internal diame- 
ter, providing a magnetic field of 3.8 T. Within the field volume are the silicon pixel and strip 
tracker, the crystal electromagnetic calorimeter and the brass/ scintillator hadron calorimeter. 
Muons are measured by gas-ionization detectors embedded in the steel return yoke in the pseu- 
dorapidity range \rj\ < 2.4, (tj = — ln[tan(0/2)], where 6 is the polar angle of the trajectory of 
a particle with respect to the direction of the counterclockwise proton beam) using three de- 
tector technologies: drift tubes (DT) (for the range \tj\ < 1.2), resistive plate chambers (RPC) 
(for \tj\ < 1.6) and cathode strip chambers (CSC) (for 0.9 < \rj\ < 2.4). The DT and RPC 
are indicated as the central "barrel" while the CSC comprises the "endcaps". A more detailed 
description of the CMS detector can be found in ||22| . 

We search the dimuon invariant mass distribution between 5.5 and 8.8 GeV (defined as "mass 
range 1") and between 11.5 and 14 GeV ("mass range 2") for a narrow resonance a, with a decay 
width ~MeV, natural in the NMSSM context. We avoid the range between 9 and 11 GeV be- 
cause the abundant contributions of the bottomonium resonances to the mass spectrum makes 
this search unfeasible. Selection criteria are applied to reduce backgrounds from QCD contin- 
uum, and we perform a mass scan in mass ranges 1 and 2 to determine a potential contribution 
from an a signal. Given the better mass resolution in the barrel part of the detector than in the 
endcaps, we also separate the mass scan into two acceptance regions, based on the dimuon rj, 
in order to improve the sensitivity. 

We analyze events collected with an online selection that requires the detection of two opposite- 
sign muons with transverse momenta pj > 3.5 GeV and additional requirements imposed at 
the high level trigger (HLT). All three muon systems, DT, CSC, and RPC, take part in the trigger 
decision. A good primary vertex is also required, as defined in Ref. [23 j. The additional HLT 
requirements include pj(}i + }i~) > 6 GeV, 5.5 < m fl + }l - < 14 GeV, and a distance of closest 
approach of the muon tracks to the beam axis compatible with that expected for prompt decays. 
A prescale factor of two was imposed on the trigger to maintain a reasonable trigger rate. 

The main backgrounds arise from QCD processes and, in the lower invariant mass range, from 
a residual tail of the Y(1S) resonance. We determine the background shape in the invariant 
mass directly from data, and use simulated events as a cross-check. Signal samples, QCD, and 
Y resonances are simulated with PYTHIA 6.4.24, Tune D6T [24J, and CTEQ6 parton distribution 
functions |25|. Tune Z2 gives compatible results. As the NMSSM is not fully implemented in 
PYTHIA, we generate the MSSM pseudoscalar A boson in the mass range of 5.5 to 14 GeV and 
require dimuon decays. These samples also contain a simulation of the effects on the number 
of primary vertices from overlapping pp interactions in the same bunch crossing. 

To select the best dimuon candidate in each event, quality criteria are applied to the tracks 
which reject misidentified muons and muons from kaon and pion decays. Muons are required 
to be within the geometrical acceptance (\t]\ < 2.4) and to be in the plateau of the trigger ef- 
ficiency, with pj > 5.5 GeV. Muon tracks are required to have at least 11 hits in the silicon 
tracker, at least one of which must be in the pixel detector, and a track fit x 2 / dof < 1.8. This 
value is chosen to maximize the signal significance with respect to the QCD continuum, which 
is extracted directly from data. 

Isolation requirements suppress misidentified leptons from jets and non-prompt muons from 
hadron decays. Muons are required to be isolated within a cone of radius AR = yj (At]) 2 + (Acp) 2 
= 0.3 around the muon direction, where <p is the azimuthal angle. The muon isolation, 7 re i, is 
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defined as the sum of the px (as measured in the silicon tracker) and transverse energy (as 
measured in the calorimeters) of all objects within this cone (excluding the muon itself), di- 
vided by the muon pj as measured by the tracker. We require 7 re i < 0.2. This requirement is 
optimized by comparing the simulated a signal with opposite-sign dimuons from data, and we 
verify that this value is appropriate for both the barrel and endcap dimuon pairs. This isolation 
requirement rejects a large fraction of the background arising from QCD production of jets. 

Dimuon candidates consist of two opposite-sign muons [26 1 with an invariant mass between 
5.5 and 14GeV. If more than one dimuon candidate is present, that with the highest x 2 proba- 
bility associated to the kinematic fit of the dimuon vertex is retained. 

The invariant mass spectrum in the search range has two main contributions: the QCD contin- 
uum and the bottomonium resonances. To characterize these shapes for use in the mass scan, 
we perform a binned maximum likelihood fit to the total invariant mass distribution. For the 
QCD continuum, we use a first-order polynomial probability density function (PDF). Each Y 
resonance is parametrized via a double Crystal Ball (CB) function p7ll28| . The resolution of 
one of the CB functions is left free in the fit but is constrained to be the same for all the three 
resonances. The resolution of the second CB function is determined from the fit of the Y(1S) 
peak, and forced to scale with the mass of the other two resonances. As the resonances overlap, 
we fit for the presence of all three Y states simultaneously using three double CB functions. 
The mean of the CB of the Y(1S) is left free in the fit to accommodate a possible bias in the 
momentum scale calibration. The number of free parameters is reduced by fixing the Y(2S) 
and Y(3S) mass differences, relative to Y(1S), to their world average values HI. 

The fits to the Y shape and continuum background are performed in the barrel and endcap 
regions separately, shown in Fig. [Tj The fitted numbers of events are given in Table [lj the 
barrel/ endcaps ratio for the Y peaks is consistent with MC predictions. Outside the Y peak 
range, corresponding to the signal search mass ranges, the data are well-described by a first- 
order polynomial. Figure [T] also shows hypothetical signals from pseudoscalar Higgs bosons 
a with masses 7 and 12GeV, and 2pb cross sections, scaled by a factor of 10 for visibility. 
We perform mass scans of the invariant mass spectra, dividing mass range 1 into 110 steps 
and mass range 2 into 100 steps of 30 MeV each, and treating the barrel and endcap spectra 
separately. At each step, we build a signal Gaussian PDF with a mean fixed to the center 
of the step and a width determined by the mass resolution, use a first-order polynomial to 
characterize the background, and perform an unbinned maximum likelihood fit to search for a 
possible contribution from the a. For each signal mass point, we determine the resolution by 
fitting the a invariant mass spectrum with two CB functions (as for the Y, the sum of two CB 
functions better describes the resolution) and the mass resolution is calculated as the weighted 
average of the widths of the two functions. The resulting dimuon invariant mass resolution 
ranges from 50 to 120 MeV (90 to 190 MeV) in the barrel (endcaps) for the mass range 5.5 to 
14 GeV. These agree well with the resolution obtained from the Y resonances in data and MC. 
We fit the resolution as a function of mass using the simulated signal samples, and use this 
to extract the values of the dimuon mass resolution for each mass bin needed in the scan to 
determine the upper limit. 

In mass range 1, we take into account the radiative tail of the Y(1S) by including its shape 
determined from the full invariant mass spectrum fit. No significant discrepancy with SM 
background predictions is observed, and we proceed to set cross section limits, as described 
below. 

The efficiency for the selection is factorized into three contributions, e = e acc x etrig x £ S ei/ where 
e acc is the kinematic acceptance for the a, e tr ig is the efficiency of the muon trigger, and e se i is the 
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Figure 1: Dimuon invariant mass distribution for the barrel (upper) and endcaps (lower) after 
the event selection. The invariant mass distributions are fitted accounting for the three Y reso- 
nances and QCD continuum. Hypothetical signals from pseudoscalar Higgs bosons a at 7 and 
12 GeV are shown. 
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Table 1: Fitted numbers of Y and continuum background events in the invariant mass range 
5.5-14 GeV. The Y contributions are summed over the three resonances. 



Contribution 


events 


(barrel) 


events (endcaps) 


Y 


93753 


±396 


95876 ± 454 


Continuum 


41210 


±320 


45792 ± 385 



efficiency of the selection applied to the dimuon candidates. We use pythia 6 to simulate the a 
signal and to determine e a cc- The trigger and selection efficiencies (e tr i g and e se j) are measured 
with J/ xp events in data using the tag-and-probe technique [ 28 1 . We perform this study in bins of 
tj and px of the probe muon. The efficiency values extracted from data are compared with those 
obtained from the simulation of prompt J/ip — > ■ The difference between the efficiency in 
data and Monte Carlo (MC) simulation is evaluated in bins of pj and rj and used as a correction 
to weight the MC events in order to accommodate possible discrepancies. These corrections are 
typically on the order of a few percent. For each dimuon candidate, the weight is the product 
of the corrections for the two muons. 

The isolation requirement efficiency that contributes to e se i cannot be measured using the ]/rp 
dataset as one of the main production mechanisms for ]/xp is through B-meson decays, resulting 
in non-isolated muons. This is not well accounted for in simulation, and would result in biased 
data/MC efficiency corrections. In order to estimate this correction, we use Z — > events 
and consider the lower pj spectrum of the probe muon. 

The total efficiency e is defined for each a mass sample as the fraction of generated signal 
events, weighted by the appropriate data/MC corrections, that satisfy all the selection require- 
ments. This ranges from 1-3.5% for the a mass range of 5.5-14 GeV, and we fit the e distribu- 
tions with second (third) order polynomial functions in the barrel (endcaps) to use in the mass 
scan. The increase in the efficiency as a function of the invariant mass is mainly due to the pj 
requirements on the muons at the HLT level. 

Several sources of systematic uncertainty affect these results, including a 2.2% uncertainty on 
the integrated luminosity p9|. The efficiency corrections are determined using the tag-and- 
probe results described above. We determine, event-by-event, the uncertainty on the total effi- 
ciency corrections by propagating the uncertainties on the single muon corrections. This total 
event efficiency uncertainty is largely independent of mass, with a maximum value of 12%. We 
apply this value as a systematic uncertainty for every bin in the scan. 

The isolation efficiency is uncertain at the 5% level, corresponding to the largest discrepancy 
between data and MC simulation in the entire relevant pj range. We evaluate the systematic 
uncertainty on the resolution of the a as the quadrature sum of the difference between the mass 
resolution of the a with a mass of 10 GeV and the resolution of the Y(2S) (which has the same 
mass) in MC simulation, and the difference between the latter and the mass resolution obtained 
for the Y(2S) from data. Additionally, the finite statistics for the determination of the mass 
resolution as a function of the dimuon mass contributes a source of uncertainty. We consider 
the mass ranges separately and include these systematic uncertainties in the calculation of the 
upper limit on the cross section times branching fraction. Overall, this adds a 11% (4%) effect 
for the barrel (endcaps). 

Systematic uncertainties on the background description include the shape uncertainty of the 
first-order polynomial fit of the background PDF. We fit the background with alternative func- 
tions (a second-order polynomial and an exponential function), generate MC pseudo-experiments 
using these functions, and fit the distributions using the first-order polynomial. The resulting 



6 



systematic uncertainties, from the distribution of the fitted parameters, is of the order of a few 
percent. 

No significant signal is observed, and we determine 95% confidence level (CL) upper limits on 
a ■ £>(pp — > a — > as a function of the dimuon mass using the CLs approach [30ti32| . A 

few steps at the edges of the mass scans, where the fitting procedure has no predictive power 
on the signal shape, are not used. Figure [2] shows the upper limit results for the two mass 
ranges including the systematic uncertainties discussed above. These limits are significant in 
the context of the NMSSM, and can be presented in terms of upper limits on | cos #a I • The larger 
the value of tan f>, the stronger is the constraint. Figure |3]presents upper limits, | cos $A| max as a 
function of m ai for tan f> = 1, 2, 3, 10, 30, 50. Our upper limits are compared to an earlier analysis 
of the BaBar Y(1S) and Y(3S) data 1331, and are superior for m ai > 7.5 GeV for tan/3 = 50, 
decreasing to m ai > 6 GeV for tan /3 = 2, and are superior for all masses at tan jS = 1. Further, 
these are the first significant limits for m a > m Y (3s)- 
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Figure 2: Upper limits at 95% CL on cr ■ B(pp — > a — > in mass range 1 (upper panel) and 

mass range 2 (lower panel) including systematic uncertainties. The dotted lines correspond to 
the expected limits, and the bands correspond to 1- and 2-a level uncertainties on the expected 
limits. 



In conclusion, we performed a search for a narrow, low mass pseudoscalar a, which is produced 
by gg — > a and decays via a — > in the mass ranges 5.5-8.8 GeV and 11.5-14 GeV, using 

a data sample corresponding to an integrated luminosity of 1.3 fb -1 collected with the CMS 
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Figure 3: Upper limits on the NMSSM parameter | cos | as a function of m ai in the two mass 
ranges. The solid curves correspond to different tan /3 values: from top to bottom, tan [5 = 1, 
tan /3 = 2, tan ft = 3, tan /3 = 10, tan f> = 30, and tan /3 = 50. For each tan f> value in mass range 
1 the second, dotted curve shows the limits from the BaBar Y analysis. There are no BaBar 
limits for tan f> = 1 in mass range 1, or for any tan f> in mass range 2. The line at |cos $A| max = 1 
is equivalent to no limit. Results from CDF are not shown as they are less stringent than the 
BaBar limits. 
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